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An improved description of four-pion production in electron-positron annihilation and in r-lepton 
decays is presented. The model amplitude is fitted to recent data from BaBar which cover a wide 
energy range and which were obtained exploiting the radiative return. Predicting r-decay distribu- 
tions from e^e~ data and comparing these predictions with ALEPH and CLEO results, the validity 
of isospin symmetry is confirmed within the present experimental errors. A good description of 
two- and three-pion sub-distributions is obtained. Special emphasis is put on the predictions for 
OJTT (— > 7r"'"7r~7r°) in e^e~ annihilation and in r decay. The model amplitude is implemented in the 
Monte Carlo generator PHOKHARA. 

PACS numbers: 13.66.Bc,12.40.Vv,13.35.Dx 
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I. INTRODUCTION 

The production of four pions in r decays and e^e~ an- 
nihilation has received considerable attention, both from 
the theoretical (iLMJs), U, a la LZIi l8|, M and the experi- 



mental side 



111 . I121 . 1131 . 1141 . 1151 . 1161 . 117| . Relating the 
cross sections and rates for the four charge combinations 
(7r+7r~27r°, 27r+27r~, 27r~7r+7r° and 7r~37r'') gives impor- 
tant hints on the validity of isospin symmetry and the 
size of the isospin breaking terms. The dependence of 
the rates and the cross sections on \fQ^, the invariant 
mass of the four-pion system, and the investigation of dif- 
ferential distributions, e.g. of the two- and/or three-pion 
masses, gives information on the resonance structure of 
the amplitude. In the low region, predictions based on 
the chiral Lagrangian can be tested which, however, must 
be complemented by resonance physics in order to prop- 
erly describe the rates in the dominant region between 
1 and 3 GeV. The e+e" cross section is, furthermore, 
important to evaluate the hadronic vacuum polarization 
which in turn is essential for the precise prediction of the 
muon anomalous magnetic moment and the running of 
the electromagnetic coupling [H, [l^ . 

From the experimental side precise r data have been 
obtained by ALEPH [ll| and by CLEO ll] collabora- 
tions, which, however, are naturally restricted to \fQ^ 
below 1.77 GeV. The e^e~ cross section has been mea- 
sured by CMD2 O, [13, [ll and SND 13] (older data 
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are far less accurate and will not used in this paper) and, 
more recently, by BaBar [l^, [l3] through the method of 
radiative return which covers energies up to 4.5 GeV. 
This method, which was proposed in @, (23,[22| allows to 
use the large luminosity at i3-factories for a measurement 
of the e+e^ cross section in the region of interest. 

From the theory side the first evaluation based on 
chiral perturbation theory has been performed by Fi- 
scher, Wagner and Wess [H and applied to r decays. 
Subsequently this ansatz was extended fs'l to include p, 
fli and /o resonances, which are clearly visible in sub- 
distributions. In addition the uj-k mode was introduced, 
again predicted from the chiral anomaly [3, [^. Later 
this ansatz, slightly modified, was implemented in the 
generator EVA |21] to simulate 47r production in the ra- 
diative return [5]. As stated above, the low region 
should be best suited for a description based on chiral La- 
grangians. Combining one-loop chiral corrections at low 
g2 with resonance enhancements at intermediate ener- 
gies, precise predictions have been obtained in [^ which 
will be discussed below. 

In view of these recent theoretical and experimental 
developments, together with need for an optimal imple- 
mentation of the 47r mode into the Monte Carlo event 
generator PHOKHARA [H [H, [H [H, [H, [B, [H, [H , 
an improved ansatz for the corresponding hadronic am- 
plitude has been developed. The ansatz is largely based 
on and [H] (concer ning the w part) with model pa- 
rameters fitted to the recent BaBar results. In order to 
accommodate the p~ signal observed in we include 
a contribution, which is modeled to mimic a SU(2) gauge 
theory with the p- meson (and its radial excitations) as 
gauge boson(s). 

Our paper is organized as follows: To facilitate the 
subsequent discussion, in Section |TT] the basic definitions 
are introduced and the (well known) isospin relations be- 
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tween the amplitudes and the rates of the four channels 
are collected. The validity of these relations is investi- 
gated in Section Unl using data from e+e~ annihilation to 
predict the corresponding, experimentally measured dis- 
tributions for T decays. The ingredients for the ansatz for 
the matrix element of the hadronic current are discussed 
in Section ITVl The comparison of this ansatz with e~^e~ 
data and the fit of its parameters is presented in Section 
IVl together with the comparison between the model and 
data for a variety of distributions. The implications of 
the model for r decays is discussed in Section IvH the 
implementation into the generator PHOKHARA and re- 
lated technical tests in Section IVIII A brief summary 
and our conclusions are given in Section IVlIII A detailed 
description of our model with the complete list of param- 
eters can be found in the Appendix. 



II. GENERAL PROPERTIES OF THE FOUR 
PION ELECTROMAGNETIC CURRENT 

General properties of the four pion electromagnetic 
current were investigated in where it was shown 
that assuming isospin symmetry just one function de- 
scribes all four matrix elements. We use the same letter 
J for the operator and its matrix element. For conve- 
nience we recall the definitions and notation introduced 
in ii. 

Ignoring the issues of isospin breaking and radiative 
corrections, the electromagnetic current can be decom- 
posed into an isospin singlet piece and a part transform- 
ing like the third component of an isospin triplet: 



jem J' 



1 



J 



1=0 



(1) 



whereas the charged current generating r decays is given 

by 



(2) 



Final states with an even number of pions are produced 
through the isospin one part only 

(7r+7r-7r?7rO|j3|0) = J^(pi,P2,P+,P") 



The function = J/^ (91, 92, 93, 94) is symmetric (anti- 
symmetric) with respect to the interchange of qi and q2 
(qs and (74). 

In it was shown that also the matrix element 
(7rj~7r^7r"'"7r''| J~ |0) can be used as an independent func- 
tion, through which the other ones can be expressed. If 
one denotes 



(vr-^^gTra^l J;|0) = F^ipi,p2,p-,P3) , (4) 
one gets relation 

(^+7r+7r-7r-|j3|0) = 

from which it is clear that the matrix element 
(TTj^vr^TTf TT^I J^|0) can be expressed by the matrix ele- 
ment (7r^7rj7r27r§| J^|0). The opposite is also true and 
we have proved it using the method developed in [6| to 
express (Tr+Tr-Tr^Tr^l J3|0) by (7r]"7r2"7r+7r°| J" |0). How- 
ever, as in both cases the obtained inverse relations are 
far from being as elegant as the ones of Eq. ([3]) and Eq. 
([5]), the result is not presented here. 

The currents defined in Eq.Q contain the complete 
information about the hadronic cross section through 



(6) 



where R{Q^) is equal to (T(e+e hadrons){Q^) / Gpoint, 
with apoint = 47ra^/(3(5^), and dl>„((3; qi, . . . , g„) de- 
notes the n body phase space with all statistical factors 
included. 

The amplitude describing the r decay into an arbitrary 
number of hadrons plus a neutrino (excluding radiative 
corrections) is given by 



Mr^ ^ V {p,) r (1 - 75) u (pr) J- , (7) 



with 



(7r+^2+^r^2-|J;^|0) = 

J^iijpt ,P2 ^Pl ^Pl) + Jt^{Pl ^P2 ^Pt ^Pl 
+ Jt,{pt-,Pl,Pt^P2) + Jt^iPl-.Pl^Pt^P2 



Jf^{P2,P3,P^,Pl) + Jf,{Pl,P3,P~,P2) + J^{pi,P2,P~,P3) 



(7rr^2-7r+^"|j;|0) = 

Jf^{P^,P2,Pl,P°) + J^{p'^,Pl,P2,P°) 



(3) 



J- 



J-(9i,...,g„) = (;i(gi),...,%„)|J-(0)|0) (8) 



and J^(0) = d'jaU at the quark level, where we restrict 
our considerations to the Cabbibo allowed vector part of 
the hadronic current. 

The differential r decay rates are given by 

fi^r — >i'-\-hadrons 



= 2 r, 



1-^) (1 + 2^) R^Q') ,(9) 
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with 

^ (QmQ- - .9p.Q') R"iQ^) (10) 

and Te = G|m^/(1927r3). Note the relative factor of 2 
between the definitions in Eq.® and Eq. fTU)) . We have 
included here also electroweak correction factor Sew (we 



use Sew = 1-0198 f?\) to account for standard elec- 
troweak corrections. 

The function is related to the spectral function de- 
fined by CLEO [11] through - +0) = 5ttV^''''° and 
through = 3vi to the vector spectral functions defined 
by ALEPH 15]. In this paper we will use normalization 
of the spectral functions chosen by ALEPH. 

The four pion spectral functions and the cross sections 
can be expressed as linear combinations of two integrals 



2n f 

A = -g^ J {qi,q2,q3,q4)J^{qi,q2,q3,q4) d$4(Q;gi, . . . ,94) 

47r f 

B = J Re (J^ (51,52,93,94) (91,93,92,94)) d$4(Q;9i, • ■ • ,94) ■ (11) 



The relations read 



R{+ - 0) ^ ^ A ; i?^( +0) = A+^B 

R{+ + - -) = A + B ; i?^ (- 0) ^ ^{A + B) . (12) 
The additional contribution to i?(-|- -I- ) — ^ + B + C of the form 

^^~Q2j (J^ (91,92,93,94) J;i(93,94,9i,92)) rf$4(Q;9i, ■•• ,94) , (13) 



vanishes for any symmetric phase space configuration. 
Eqs[T2]correspond to the familiar relations between r de- 
cay rates and e+e~ annihilation cross sections: 

i?^(- 0) = \r{+ + - -) 

i?"(- - +0) = i i?(+ + - -) 

+ R{+ - 0) . (14) 

III. ISOSPIN SYMMETRY - EXPERIMENTAL 
SITUATION 

In this section we would like to address the question, if 
present experiments require inclusion of isospin violating 
effects in the model. Combining the results from BaBar 
[1^ on a{e'^e^ 2it'^2'k~) with their preliminary results 
on (T(e+e~ 27r°7r+7r~) |l3| and using Eqs.(ITl|), one 
obtains predictions for the r spectral functions. These 
can be compared with ALEPH [ll| and CLEO [ul data 
(compare also 17]). As shown in Fig. [1] and Fig. [21 
T- and e^e~- data are in good agreement within the er- 
rors, even if one observes systematical shifts. However, 



these shifts are well within the 5% systematic error of 
CLEO and the 6% (27r~7r+7r") and 10% (tt-Stt") errors 
for ALEPH spectral functions (ALEPH does not give sep- 
arately the systematic error) as well as the 5% to 12% 
systematic error for BaBar cr(e+e" — > 27r+27r^). For the 
preliminary BaBar data [13] on (T(e+e^ 27r"7r+7r~) a 
10% systematic error is assumed. Truly isospin breaking 
effects are expected to occur at thepercent level due to 
the TT^ — tt" mass difference alone [5|. 

From the cross section a{e^e'^ 27r°7r+7r~) and the 
relative contributions of the wtt final state as given in [l3] 
(since the errors were not specified there, we attribute 
20% error to the spectrum) one can infer the ld contri- 
bution to cr(e+e^ 27r'^7r+7r^). Based on this result 
one can predict the omega part of the — > i^27r^7r+7r*' 
spectral function and compare it with the CLEO result. 
Satisfactory agreement is observed in Fig. [3] 

From the comparisons of experimental data we con- 
clude that no isospin symmetry violation is observed 
within the present accuracy. Thus the model we propose 
to describe the data is based on isospin symmetry. How- 
ever, effects from the pion mass difference in the phase 
space are included. 
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FIG. 1: (color online). The spectral function of the r~ — > 
u3'k'^tt~ decay mode. ALEPH data [l^ versus predictions 
from BaBar data [l^. [iTj and the model predictions. 




0^ (GeV) 



FIG. 2: (color online). The spectral function of the r — > 
u2-k'tt+7v° decay mode. ALEPH [lland CLEO [H data 
versus predictions from BaBar data [id. [l7| and the model 
predictions. 



IV. THE MODEL OF THE FOUR PION 
ELECTROMAGNETIC CURRENT 

There are many motivations why the model adopted 
in should be updated. First of all new and more ac- 
curate data are available. The CLEO data on tau de- 
cays [ni, which were not used in the tau spectral 
functions from ALEPH fis'l and the measurement of the 
cross section of the reaction e+e~ 27r+27r~ via the 
radiative return method by BaBar 16] provide us with 
the opportunity for a substantial improvement of the 
model implemented in the event generator PHOKHARA 
m m, Hi, m, IMJp, m, ^ . The omega part of the cur- 
rent, which in [3|, |5| was implemented without structure, 
is now known much better from phenomenological stud- 
ies [2^. The new preliminary data from BaBar [17] on 




FIG. 3: (color online). The omega part of the spectral func- 
tion of the T~ — > //27r~7r"'"7r° decay mode. CLEO lllll data 
versus predictions based on preliminary BaBar data [l7| ] and 
the model predictions. 



than implemented in . All this was taken into account 
in constructing the model presented in this paper. 




the reaction e+e 



27r'^7r+7r also show richer structure 



FIG. 4: Diagrams contributing to the hadronic current in 

The amplitude used in Q is schematically depicted 
in Fig. m In the contributions from the first two di- 
agrams, which proceed through the intermediate reso- 
nances p — > fliTT and p ^ fop respectively (where p 
stands for p(770) and its radial excitations), only the 
parameters of the current were adopted to the improved 
data and a new p resonance (/9(2040)) was added (neces- 
sary to fit the BaBar fld\ data). The contribution from 
uj, where previously the substructure of the omega decay 
was not taken into account, is now modeled using infor- 
mation from (28j . Schematically, the new cj amplitude is 
depicted in FigO 

BaBar has, furthermore, observed [TZj a strong p^p^ 
contribution. Thus a new part containing the p pp 
contribution has been added, treating the p particles like 
SU{2) gauge bosons. The contributions to the amplitude 
are depicted in Fig. [G] For more general frameworks, 
where such terms are present see [30| (and references 
therein). The SU(2) symmetric Lagrangian describing 
p— pair production reads 
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FIG. 5: The new contributions from the omega part of the 
current. 





FIG. 6: The new contributions from p mesons. 



c„ = 



4 
1 2 
2™^ 



fiiy 



■ F 



1 2^ 
2^pP^ ■ P 



(15) 



where 




and 



9Pa X Pv 




D„ 



df.^ +g[ Pn^ 



(17) 



The only free parameter, the couphng constant g {g = 
gpTTTr), can be extracted from p tttt decay. However, 
as it stands, the model leads to a wrong high energy be- 
havior of the cross section, falling less rapidly then the 
data. This problem can be cured by adding p' contribu- 
tions and allowing for trilinear couplings between p and 
p' . It was also necessary to relax the fixed coupling g to 
fit the data. The detailed description can be found in the 
Appendix. The model can be further refined, when more 
experimental information is available. 

The behavior of the four pion amplitude in the low 
region has also been studied Q in the framework of 
chiral resonance theory, including terms up to 0{p'^) [3l| . 
The implementation of resonances and their parameters 
differs from the choice in this paper. The results of the 
two models are compared to the data in Figs. [TUlandfTT] 



FIT OF THE CURRENT PARAMETERS TO 
THE EXPERIMENTAL DATA 



To separate the well measured to contribution, from the 
rest, we fitted the parameters of the model to the lo part 
of the cross section of the reaction e+e~ 27r"7r+7r^ ex- 
tracted from preliminary BaBar data [r?|. Furthermore, 
we fitted the model parameters to the cross sections of 



measured by BaBar [l^ [T3| • 




Q^{GeV) 



FIG. 7: (color online). Fit to the u— part of the cross section 
for e^e~ — > 2-k^-k'^'k~ (20% systematical error for the 
preliminary BaBar data was assumed). 

The results are shown in Figs. [71 [S] and [H and in Table 
im The fit is quite good, with xV^d.o./ = 275/287. How- 
ever, one has to remember that only the cross sections 
were fitted and all sub-distributions are to large extent 
determined by the model assumptions. The constants 
and f3f (with i = 1,2,3) characterize the rela- 
tive importance of the radial p excitations (compared to 
the one of the ground state, /3o = 1) in the amplitudes de- 
picted in Figs, m and[5](see also Eq. (jA.ll[) ). Large values 
of I3(° are the consequence of the small /o (1370) — po — Po 
coupling compared to higher p radial excitations, which 
are indeed dominated by pi. 

It is interesting to see how the model compares to pre- 
dictions based on the chiral Lagrangian [6] in the low 
region, where this ansatz is expected to be applicable. 
In Fig. [TO] (Fig. \TT^ this comparison is shown for the 
charged (neutral) mode, together with data from BaBar 
[11 [i3, CMD2 [H m El and SND [l|. Since our 
model parameters were fitted to that cross section the 
thick dotted curve is not a prediction, apart of the ^/Q^ 
region below 0.8 GeV, from where the contribution to 
of the fit is negligible due to the low accuracy of the data. 

The sub-distributions can be qualitatively compared 
(Fig. [T^ with plots presented by BaBar These 
were not used in the fit and thus can be considered as 
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FIG. 8: (color online). Fit to the data for a{e^e~ — + 
2tv''tv'^'jt~), taken from [ITI ] (10% systematical error was added 
to the statistical error). For comparison also CMD2 [T^] and 
SND 13] data, which are consistent with BaBar data, are 
shown (without their 10%-20% error-bars). Contributions 
from p part of the current (Eq.(|Aj8|) to the cross section 
(see text for definition) are also shown. 
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FIG. 9: (color online). Fit to the data for a{e^e~ —> 
27r"'"27r~), taken from [la ]. For comparison also CMD2 
[lol . Il4| and SND [l3] data, which are consistent with 
BaBar data, are shown (without their 7%-20% error-bars). 



predictions. Integrals for both experimental and theo- 
retical plots are equal by construction. Further refine- 
ments of the model will be possible when the data on 
sub-distributions will become available. 

The contributions from two p mesons in the final state 
are shown as dashed line in Fig. \8\ They were extracted 

7r° invariant masses 



FIG. 10: (color online). Comparison of our model with the 
chiral Lagrangian predictions (solid and dotted lines) in 
low energy region and experimental data [13. [l^. . 
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FIG. 11: (color online). Comparison of our model with the 
chiral Lagrangian Q predictions (solid and dotted lines) in 
low energy region and experimental data [lol, W\ . 



affected by background from the other amplitudes and 
thus do not correspond exactly to the contributions from 
p part of the current (Fig. [S] and Eq. (|A.8[) ). hence the 
separation is not as clean as for the w case. The model 
prediction is smaller than the BaBar result [iTj . 

Selected two and three pion invariant mass sub- 
distributions for the reaction e+e^ 27r'^7r+7r~7(7) are 
shown in Fig. 1131 The contributions from various reso- 
nances included in the model are clearly visible. Com- 
parisons of the predictions will be possible, when the final 
BaBar results are published. 
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FIG. 12: (color online). Two and three pion invariant mass distributions for five difi^erent ranges of 2-k'^2tv~ invariant mass. 
The BaBar data points (filled circles), given as events/bin, are superimposed on plots obtained by PHOKHARA (filled squares) 
(see text for details). 
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VI. MODEL PREDICTIONS FOR r DECAYS 

One can confront the model with the data ^33] for the 
partial t decay rates to four pion final states. The results 
are collected in Table HI The theoretical error is obtained 
from the errors of the model parameters extracted in the 
fit. Within the quoted errors, the predictions are in good 
agreement with the data even if one observes sizable dif- 
ference between the data for Br{T^ Vt-Itt^tt^'k^) and 
the prediction via the isospin relations. At present the 
results are still consistent within the conservatively esti- 
mated error, which is dominated by the one of the pre- 
liminary BaBar result for a{e^e~ — > With 
an expected error of about 5%, the final BaBar result 
will further push the accuracy of the isospin symmetry 
tests. 



The model of 47r hadronic current proposed in this pa- 
per was fitted to BaBar data and relies on isospin sym- 
metry. Thus its predictions for the r spectral functions 
follows the predictions from BaBar data based on the 
isospin symmetry assumption (presented in Section 
apart from few percent phase space effects coming from 
the TT^ — TT*^ mass difference. The model predictions are 
also shown in Fig. [1] [2] and [31 The central curve repre- 
sents the model predictions, the upper and lower curves 
are the error estimates based on errors coming from the 
fitted parameters of the model. 

Two, three and four pion invariant mass distributions 
obtained within our model are compared with CLEO 
data (available only as plots) in Fig. [T4l Although pre- 
dictions and the data differ as far as the detailed descrip- 
tion is concerned, good qualitative agreement is observed. 





Br(r > Vr2n- Ti+ n^) 


Br(r > j^rTT a;(7r vr^'^Tr'^)) 


Br(r > v-r-K^Z-K^) 


PDG [32] 


(4.46± 0.06)% 


(1.77± 0.1)% 


(1.04± 0.08)% 


model 


(4.12 ± 0.21)% 


(1.60 ± 0.13)% 


(1.06 ± 0.09)% 


BaBar (CVC) 


(3.98 ± 0.30)% 


(1.57 ± 0.31)% 


(1.02± 0.05)% 



TABLE I: Branching ratios of t decay modes. Results of our model are compared to experimental data [33| and predictions 
based on BaBar data [3, [l^l and isospin symmetry 



VII. IMPLEMENTATION INTO PHOKHARA 
AND TESTS OF THE MONTE CARLO 
GENERATOR 

The model for the hadronic current was im- 
plemented into the PHOKHARA event gen- 
erator (version 7.0). It will be available at 
|http : //if ic . uv . es~rodr igo/phokhara/ together 
with the implementation of the J/ijj and ijj{2S) contri- 
butions to 2-body hadronic final states (in preparation). 
Only the current was coded in the form described 
in the Appendix; the charged mode is obtained via the 
relation Eq. Neither the uj- part of the current 

nor the double p resonance diagrams (left in Fig. [6]) 
contribute to that part. Only a priori weights, used in 
the multi-channel Monte Carlo generation, were changed 
as compared to previous versions (2^. Nevertheless, 
tests checking the implementation were performed to 
assure a proper technical precision of the code. The NLO 
version of the code was checked for the configurations 
without any cuts against analytic NLO results of [s^ 
(see also [2j]), separately for one- and two- photon 
contributions. The separation w = = 10^^ between 

soft (integrated analytically) and hard (generated) parts 
was used in this test. The precision of the tests, limited 
by the Monte Carlo statistics, is significantly below one 
per mill. As the analytic formula contains as a factor 
the cross section of the process without photon emission 



{a{e'^e~ — > 47r)), which is not known analytically, it 
was obtained be means of Monte Carlo integration by 
a dedicated program. In that program, in contrast to 
PHOKHARA, flat phase space generation was used to 
avoid any errors due to the change of variables. The 
independence of the results of the generation on the 
separation into soft and hard parts was also tested with 
similar precision. 



VIII. CONCLUSIONS 

Four-pion production in e+e~ annihilation and in t 
decays is characterized by a multitude of resonant sub- 
channels. This makes it difficult to construct an ampli- 
tude which is based on first principles of QCD only. In 
this paper we have constructed a model amplitude which 
incorporates a limited set of channels, namely ai7r, pfo, 
pp and LOTT and which is approximately consistent with 
chiral predictions for small Q^. A number of parameters 
which characterize the relative importance of the various 
couplings and of the radial excitations of the p meson is 
fitted to the cross sections for 27r+27r^ and 7r+7r^27r° as 
measured by the BaBar collaboration. The model pre- 
dictions for the two- and three-pion mass distributions 
which are not fitted separately, are consistent with the 
data both from e"'"e~ annihilation and from r decays. 
Furthermore, we find that the present data are, within 
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FIG. 13: (color online). Predicted by the model two and three pion invariant mass distributions for five different ranges of 
27r°7r"''7r~ invariant mass (selected set). 
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FIG. 14: (color online). Two and three pion invariant mass 
sub-distributions for different pion charge combinations. The 
CLEO data points (filled circles), given as events/bin, are 
superimposed on plots obtained within studied model (filled 
squares). See text for details. 



their 5 - 10% systematic error, consistent with the re- 
lations derived from isospin invariance, and which are 
intrinsic for our model. The amplitude is incorporated 
into the Monte Carlo generator PHOKHARA, simulating 
47r production through the radiative return. 
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APPENDIX: THE CURRENT 

In this appendix we give a complete definition of the 
hadronic current used in this paper. We take as a ba- 
sic building block the electromagnetic current J^"* = 
^ (using J^"° — 0), which is of direct importance 
for implementation in PHOKHARA and define r** = 
--^(7r+7r~7r°7r2 1 J^|0). Other channels can be obtained 
through Eq.dS]). 

can be decomposed into the following four parts: 



ai 



p" 



(A.l) 



We denote the four pion momenta by 
92(7r°), '73(7r^) and q4{TT^) and use the proper 
pion masses m^± and m^o wherever appropriate. Thus 
the current possesses isospin symmetry, broken only by 



kinematic effects. The general structure is largely based 
on [5|. Contribution from the part containing an ai 
exchange reads: 

Tai ('71,92,93,94) = 

f Qi (93 , 92 , 91 , 94) + f Qi (93 , 91 , 92 , 94 ) 
- f{j\((j4,g2,9i,93) - f(^i(94,9i,92,93)- (A.2) 

The function F^^ is of the form 

f^^,(9i,92,93,94) = c„, Fp(Q^J-') BWa,{{Q-qif) 

192(93 - 94) 



;Sp((93 + 94)' 
-Q 



(93 - 94)^ + 9'i 



(g-9i)2 



^ f Qjqs - 94) (Q9i)(92(93 - 94)) 



QHQ~qiy 



(A.3) 



and corresponds to the configuration p{Q) — > 7r(gi) ai(— > 
p 7r((72)) with p T^iqa) T^iqi)- The other three terms 
are enforced by Bose symmetry {qi 52) and charge 
conjugation. 

The contribution from p /o(7''(9i) ''■(92))p(^ 

•^(93) T^iqi)) reads 



r;„(9i,92,93,94) = Cf, Fp (q2,/3-^«) ((93 + 94)') 
BWf, {{q, + q,f) \{q, q.Y ~ Q'^^"' " 



(A.4) 



The contribution coming from the anomalous part of 
the current (containing uj exchange) reads 

r^i(9i,92,93,94) = fj;;(9i,92,93,94) + fi;(92,9i,93,94) 

(A.5) 

with 

fi;(9i,92,93,94) = 2 c„ g^T,p gp^^ Fp{Q^,(3'^) 
X BW^iiQ - qif) Hp{{q2 + q^f, (92 + 94)', (93 + 94)') 
X [92 ((9i94)(93<9) - (9i93)(94(5)) 
+ 93 ((9i92)(94Q) - (9i94)(92Q)) 
+ 94^((9i93)(92Q) - (9i92)(93Q))]. (A.6) 



where 



42.3 GeV-^ 



.9p7I 



5.997. 



(A.7) 



The first term in Eq. ljA.Sp corresponds to the configura- 
tion p — > 7r((7i)cj(-^ 7i'(92)7i'(93)7r(94)) the second follows 
from the Bose symmetry. 

The structure of the omega decay was taken from [2I] . 
Other possible contributions to the Stt part of the current 
coming from (/)(1020) or higher radial uj excitations are 
not seen in the data and thus were not included into the 
model. 
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part of the current is of the following form: 



rp ((71,92,93, 94) 



,92,93,94) + G^(94,9i,92,93)) - (3 4) 



1 



} 



(A.8) 



where 

(91, 92, 93, 94) = q^BWP^-P^ ((91 + 93)') 



BWP°-P' ((92 + 94)')(0 + 293)(92 - 94) + 2 (A.9) 



and 



BW^{p\mp,Tp)lml^~BW:i{p',mp„Tp,)/ml^ . 

(A.IO) 

For the p — 7* coupHng we use gp^ — 0.1212 GeV^. 
The double resonant terms disappear in the (-000) and 
(-I-H — ) channels, the single resonant contribution, how- 
ever remains. 

For completeness we list all propagators required for 
the current. A new contribution was included in 
Fp iff') only. The different p propagators Tp, Fp and 
Bp are used in the current due to the fact that the p may 
couple in a different way to different resonances and the 
propagators themselves contain indirectly some informa- 
tion about the couplings. 



1 



+ /3i + /32 + /33 



[BWi{Q\mp,Vp) 



PiBWsiQ^ fhp, , f p J + P2BW3 (Q2 , rhp, , fp, J 



+[33BW3{Q\fhp,,fp,)] , 

where f3 = {l3i, P2, Ps) and 
W3(Q^mp,^p) = 



(A.ll) 



iTpTTlp 



(A.12) 



Q^-4m^ 



— 4m^ 



Only the masses fhp- and the widths Fp^ of pi , p2 and ps 

that appear in Fp (^Q^,P^ were fitted to the data. The 

results are listed in Table |TT1 

For the masses and the widths of particles in all other 
parts of the current we use their PDG values: 

TOp = 0.7755 GeV, Fp = 0.1494 GeV. 



rup, = 1.459 GeV, Fp, = 0.4 GeV, 
= 1.72 GeV, Fp^ = 0.25 GeV. 



(A.13) 



nip. 


1.437(2) 




1.738(12) 




2.12(2) 




0.520(2) 


r 


0.450(9) 


Fp3 


0.30(2) 




-0.066(3) 




-0.021(1) 




-0.0043(5) 




7(6) -lO-* 


Pi" 


-2.5(5.0) -10^ 




1.9(1.6) -10^ 




-0.33(8) 




0.012(3) 




-0.0053(7) 




-225(3) 




64(3) 




-1.47(4) 


Cp 


-2.46(3) 




275 


nd.o.f 


287 



TABLE 11: Values of the couplings masses and widths ob- 
tained in the fit. Masses and widths in GeV; couplings f3^ , 
{j = ai, /o, w and i = 1, 2, 3) as well dimensionless; 
couplings and c/^ in GeV~^; coupling in in GeV~^. 



Bp{Q^) - [W3(Q^mp,Fp) 

+ /3W3(Q',mp,,FpJ]/(l + /3),(A.14) 

with (3 = -0.145 . (A. 15) 

The oi propagator reads: 

,^2 



with d, [33 



(A.16) 



5(Q2) ^ 1.623 Q2 _|. 10.38 - 



9.32 0.65 



Q2 (Q2)2 

for > {mai + m^)^, 

g{Q^) = 4.1 {Q^~9mlf[l~i.?,{Q^~%ml) 
+ 5.8 (Q2_ 9^2)2 

for Q2 < [ma, +m^f (A.17) 

(g2 in GeV2) and 

ma, = 1.23 GeV, F„, = 0.2 GeV. (A.18) 



Tp{Q^) = [BW3iQ\mp,rp)+PiBW3{Q\mp„rpJ 
+ p2BW3iQ^, mp,,Tp,)] /{! + p,+ p2){A.19) 



where 

Pi = 0.08, 02 = -0.0075 . 
The /o meson propagator chosen to be 



(A.20) 



m 



fa 



(A.21) 



12 



where 



= 1.35 GeV, r/o=0.2GeV 



(A.22) 



HpiQlQlQl) = BW3{Qlmp,Tp) + BW^iQlmp.Tp) 
+ BW3{Qlmp,rp). (A.23) 



is the a; meson propagator with 



(A.24) 



0.78265 GeV, = 0.00849 GeV. (A.25) 

To obtain the correct chiral hmit [6] of the current 
Eg. ljA.ip the following relations should hold 



4 1 
"3/1 

4cp gl^^ gpi 



CA.26) 



where /^r = 0.0924 GeV. Comparing the fitted value 

c^J = -225(3) GeV"^ (A.27) 
with its proper chiral limit 



-156 GeV" 



(A.28) 



and also the second relation from Eq. (jA.26P as obtained 
in the fit 



2^1 



PO 



Pi 



-273(5) GeV"^ 
-307(4) GeV"2 



(A.29) 



it is clear that they hold only approximately. This reflects 
the fact that the fit was performed in a Q^- region, where 
the chiral Lagrangian leads, at best, to an approximate 
treatment, and that furthermore higher order terms are 
present. 
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